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Abstract 

Non-stationary column accretion onto a surface of a magnetized neutron star 
is studied with a numerical code based on modified first-order Godunov method 
with splitting. Formation and evolution of shocks in the column is modeled for 
accretion rates ranging from 10^^ g s~^ to 10^^ g s~^ and surface magnetic fields 
ranging from 5-10^^ G to 10^^ G. Non-stationary solutions with plasma deceleration 
at collisionless oscillating shocks are found. The kinetic energy of the accreting flow 
efficiently transforms into a cyclotron radiation field. Collisionless stopping of the 
flow allows a substantial part of accreting CNO nuclei to avoid spallation and reach 
the surface. The nuclei survival fraction depends on the surface magnetic field, 
being higher at lower magnetic fields. 
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1 Introduction 



Accretion onto compact objects has been being considered an efficient source of hard emis- 
sion for about forty years aheady (see the pioneer papers of Zel'dovich (1964), Salpeter 
(1964), Pringle and Rees (1972), and Shakura and Sunyaev (1973)). 

A number of successful analj^ical, semi-analytical, and numerical models of accretion 
onto neutron stars (NSs) and black holes has been developed by now [e.g. see the book 
of Frank, King, and Raine (2002) and references therein], but several important problems 
already formulated in the very earliest papers have no clear answers as yet. One of 
the fundamental points is whether the matter accretes as a gas of interacting particles 
(" hydrodynamic regime") or as a number of separate non-interacting particles ("freefall 
regime"). Zel'dovich (1967) and Zel'dovich and Shakura (1969) considered both cases and 
showed that the spectrum emitted near the surface of a NS critically depends upon the 
particular regime under consideration. Up to now accretion models exploit only indirect 
qualitative arguments to give prove to a particular regime in the hope that a realistic 
model supported by observational evidences may a posteriori justify their choice at least 
for a certain range of global parametres of an accretion system. Another important point 
is whether and when collisionless shocks (CSs) appear and evolve in an accreting flow to 
decelerate the matter on its way to the surface. 

Bisnovatyi-Kogan and Pridman (1970) pointed out that a CS can appear in a flow 
accreting onto a NS if the star possesses a dipole magnetic held B ~ 10^ G. 

A shock that decelerates the accreting flow in a binary system plays a key role in 
the model of Davidson and Ostriker (1973). Shapiro and Salpeter (1975), Basko and 
Sunyaev (1976), Langer and Rappoport (1982), and Braun and Yahel (1984) considered 
accretion onto a NS under various assumptions on the star's magnetic fleld strength and 
found stationary solutions for similar sets of hydrodynamic equations used to describe the 
accretion flow. The models of Shapiro and Salpeter (1975) and Langer and Rappoport 
(1982) postulated the existence of a stationary collisionless shock at a certain height 
over the star's surface, the height being a model parameter. Basko and Sunyaev (1976) 
constructed an accretion model with a radiative shock in the star's atmosphere. Braun 
and Yahel (1984) claimed that a stationary collisionless shock can exist over the surface 
of a magnetized NS only when the accretion rate is low enough (namely, when the rate 
does not exceed a few percent of the Eddington value). 

Detailed two-dimensional non-stationary radiation- dominated super- Eddington models 
of accretion onto a magnetized NS have been developed by Arons, Klein, and collaborators 
(see Klein and Arons (1989), Klein et al. (1996) and references therein). The authors 
of the models showed that if the accretion rate is much higher than the Eddington one, 
non- stationary radiation- dominated shocks come out and evolve in the accretion column. 
An essential feature of these models is that the existence of a shock in the column is not 
postulated a priori, being a result of evolution of the accretion system. 

In this paper a numerical model of suh-Eddington one-dimensional non-stationary hy- 
drodynamic (in the above-mentioned sense) accretion onto a magnetized NS is presented. 
The main feature of the model is a Godunov-type numerical approach that allows to 
describe discontinuous flows, and in particular, dynamics and evolution of shocks. 

The model consistently takes into account kinetics of electron-ion flows in a strong 
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external magnetic field. The field in the accretion column is fixed, but it plays a key role 
in various processes of interaction of matter and radiation that govern the flow evolution. 

In section |21 a hydrodynamic model of an accretion flow is presented. Based on the 
modeled flow profiles, destruction probabilities for accreting CNO nuclei are calculated in 
sectional The obtained results are summarized in section |3] 



2 Accretion Flow Hydrodynamics 
2.1 Formulation of the Problem 

Evolution of an accretion flow in a magnetic column above a polar cap of a magnetized 
NS at distances upto several NS radii from the surface is considered. 

The accretion flow is assumed to be "hydrodynamic" as magnetohydrodynamic insta- 
bilities grow very fast under typical conditions of an accretion column. As the initial flow 
is highly anisotropic, multistream instabilities can have growth rates comparable with 
the ion plasma frequency oOpi ~ 1.2 ■ lO^^nJg^ s~^, where Uig is the ion number density 
measured in 10^^ cm~^. Such growth rates are typical for isotropisation processes in CSs. 
Note that the cyclotron frequencies are much higher at considered magnetic fields. A 
typical propagation velocity of MHD-waves is V^jf^^^ ~ c(l — a), where a=10^^^- nig- 
^12 ^ 1- This scaling formula can be derived from general relations mentioned, e.g., in 
the book of Libermann and Velikovich (1986). Microscopic modeling of instabilities at 
NS magnetic fields is not yet feasible, and thus it is assumed that instabilities grow with 
increments of order Upi. 

It is assumed that electrons and ions move in the accretion flow with the same mean 
(flow) velocity but have different temperatures. 

The NS is assumed to possess a constant dipole magnetic field on the considered 
timescales. The accretion column configuration is shown in FigO 



2.2 Main Parametres and Equations of the Model 



The main parametres of the model are the mass (M^ ) and the radius (r* ) of a neutron 
star, as well as the magnetic field strength at its magnetic pole (B* ) and the accretion 
rate per unit area of the base of the accretion column (M/Aq). 

A basic system of hydrodynamic equations that describe flow evolution can be written 

as 
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where p = pe + Pi, P = Pi+ Vei ^ and denote sources of momentum and energy, while 
s = i,e denotes the sort of particles. 



The system should be completed with equations of state for each sort of particles. The 
equation of state for an ideal gas is used: Eg = Ps/[Ps{ls — !)]• For Xg = ksTg/ (m<jC^) ^ 1 
the adiabatic index 7^ may be written as 7^ ~ 7os(l — Xi) (de Groot et al. (1980)), where 
7oj = 5/3 is a usual non-relativistic value for particles with three degrees of freedom, while 
7oe = 3, as in strong fields under consideration the electrons are quasi-one- dimensional. 

The momentum distribution of the electrons is one-dimensional, because the typical 
relaxation time of electron Landau levels is about 10^^^Bi2 where Bi2=B/10^^ G (see 
e.g. Bussard 1980); this time is the shortest one in the system after the cyclotron time. 
The radiative decay time of the ion Landau levels is about 5-10~^Bj~2^ s that is much 
longer than the time of collisionless relaxation of the ions oc cj"/. Thus, the ions occupy 
highly excited Landau levels and may be considered three-dimensional and treated quasi- 
classically. 

One- dimensional motion of accreting plasma along dipole lines of the star's strong 
magnetic field is considered. Within such a geometry the system Q-© may be rewritten 
as 



„3^P , 9 ,3 



^'^ + ^(^» = (4) 



dt dr 



■7-[r {p + )] = r^J^ + 3ry (5) 



3^ 
r — 

dt 



Ps{Es + Y 



d ti^ 
+ —{r^lpMEs + ^)+ Psu]) = r^Qs. (6) 
or 2 



The system (H))-® should be completed by initial and boundary conditions. An 
accretion column filled with cold free-falling gas is considered as an initial condition for 
the simulations. The boundary condition in the upper part of the column is an inflow of a 
cold supersonic stream, while the boundary condition at the star's surface is the absence 
of any flows into the star. 

2.3 Physical Processes in an Accreting Flow 

In this section the processes contributing to JF and Qs in the system (jH)-® are described. 

Since a single- velocity flow is considered, the forces acting on ions and electrons are 
all summed in a single force term: 



T^-* 1 -'- nrn.i: ^ atmi ^ ^ arnii ^ nonres J- 7 



res } 



where F^^.^^ and Fatm denote gravity and viscous friction force (that is substantial only 
in the thin atmosphere of the star), while F^^^^, Fnonresi and F^es denote gravity, non- 
resonant, and resonant radiative pressure, respectively. 
The force of gravitation acting on a unit volume is 

Fgrav = + = (^"^i + Zum^)^^, (7) 

where n - is the ion number density. 
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To calculate the non-resonant radiative pressure the following formula from the book 
of Zheleznyakov (1996) is adopted: 



Up 



1 + tt 

where rig is the electron number density, (Tst is the Stefan-Boltzmann constant, is the 
local temperature of radiation field, and tt is the non-resonant optical depth. 



An equation of radiation transfer in the cyclotron line is numerically integrated to 
calculate the resonant radiative pressure force F^es = dUphot/dr and obtain Uphot - the 
energy density of the photon field. 

As the accretion column is optically thick in the cyclotron line, the transfer equation 
may be written as a diffusion equation: 
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where Sphot denotes the sources of cyclotron photons, and Jphot = |u Uphot — f^^Uphot is 
the diffusive flux of photons in the cyclotron line. As diffusion of cyclotron photons across 
the magnetic filed lines is substantially hampered (see e.g. Arons, Klein, and Lea 1987), 
only the parallel component of © is considered, that may be written as 
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where k,\\ is the diffusion coefficient parallel to the magnetic field. The equation is inte- 
grated by a shooting method with the following boundary conditions: cyclotron photons 
freely leave the upper border of the column, while their number density on the star's 
surface is determined by a blackbody spectrum of temperature Tg//. The surface tem- 
perature Teff varies with the part of flow energy released at the surface. 

To calculate the friction force acting on the flow in the star's atmosphere a standard 
expression for Coulomb stopping in a dense environment is used: 
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where A is a Coulomb logarithm, Ua is the electron number density in the atmosphere, 
and u is the flow velocity. A similar expression was used in the paper of Bildsten et al 
(1992) devoted to coUisional destruction of CNO nuclei in a NS atmosphere. 



The forces acting on the flow perform a certain work on it, and this work is effectively 
redistributed between the ions and the electrons. 

Let an external force act on the ions and an external force F^ act on the electrons, 
then it can derived from the local electroneutrality (if the frequencies of variations of the 
external forces are much lower that the plasma frequencies) that the fluxes of both types 
of particles are equal, i.e. 



where E is the ambipolar electric field, and is the effective electron-ion relaxation 
frequency. Since — Zrii, 

eE = F""— F' ^ 



where ^ = me/rrii, and the resulting effective force acting of the ions is 

F:^ = iF^ + F^Z)^, 
while the resulting effective force acting on the electrons is 



The temperature of the ions changes in the following processes: 

(i) small-angle e-i Coulomb scatterings ("heat exchange"): Hjg 

(ii) coUisional excitation of electron Landau levels: Qcyc 

(iii) collisional Coulomb relaxation in the star's atmosphere: Qreiax 

(iv) the work of the effective force: F^qU 

The model also accounts for the cyclotron cooling of the ions. This effect is noticeable 
if the star's dipole magnetic field exceeds 5-10^^ G and the column is transparent for the 
photons of the ion cyclotron fundamental line. 

The temperature of the electrons changes in the following processes: 

(i) small-angle e-i Coulomb scatterings ("heat exchange"): Hgi 

(ii) Bremsstrahlung cooling in c-i and e-e collisions: Bvei Bvee 

(iii) collisional excitation of Landau levels in e-i and e-e collisions: CyCei CyCee 

(iv) Compton processes: Qcompt 

(v) the work of the effective force: F^^u 

For a detailed description of the terms above, see the Appendix. 
2.4 The Numerical Approach 

A multicomponent accretion flow onto the surface of a NS can have discontinuities, and in 
particular, shocks. That is why the chosen numerical approach is based on a well-known 
Godunov method (see, e.g. Godunov 1959; 1964). 



The traditional form of Godunov method is only applicable to one-component systems^ 
without source terms. To account for the source terms in (0])-® a numerical scheme based 
on the work of LeVeque (1997) is exploited. The system (jlj)-® is split into two parts: 
one of them consists of fiux-conserving terms and is integrated with a modified Godunov 
method, while the other describes the sources of energy and momentum and is integrated 
with a Gear method. 

Such a way of solving the system is due to its complex structure with two sorts of par- 
ticles non-linearly interacting with each other, with external magnetic and gravitational 
fields as well as with the radiation field. 

To make the system dimensionless JF is multiplied by Cp = while Qs is 



multiplied by Cq = -y^, where t*, u^, and are scales of time, velocity, and density, 
respectively. 

Integration of equations (^-(jHl) in the accretion column is performed in a combined 
way that allows to generalize the standard Godunov method for systems with energy and 
momentum exchange between the components (i.e. account for source term in (jl))-®)- 
The accretion column is represented by a number of spatial cells. The system is integrated 
from an initial state at the moment t = to a current state at a moment t in a number 
of time steps At each of them consisting of the following stages: 

(i) integration of fiux-conserving terms without source terms: 



Mpu) d 3 2n 
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dt 



At this stage the system is simultaneously integrated in all the cells, 
(ii) integration of source terms in individual cells: 
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where qi are cell-averaged quantities in the cell number i. 

At each of the two stages the system is integrated over the same time step determined 
by the Courant condition at the first stage. 

The system (fT3j) - (fT^ is stiff and thus it is integrated with a standard LSODE routine 
(Hindmarsh 1983) that implements the Gear-B method (see, e.g. Gear 1971). It should 

^For a review of modern multi-component methods see, e.g., Zabrodin and Prokopov (1988). 
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be noted that a number of modern methods for integration of stiff systems based on the 
modified Buhrsch-Stoer method (see Press et ah 1993) require an exphcit Jacobian of the 
integrated system to be supphed. As the system (fT3j) -(fT ^ has a very complex right-hand 
part with non-analytical terms, the Bulirsch-Stoer methods appear too complicated to 
implement . 

The system (I12p - (jl4p is integrated with a "capacity-differencing" modification of the 
standard first-order Godunov method suggested by LeVeque (1997). 

The essence of "capacity-differencing" is the following. If a conservation law for a 
physical quantity q(x,t) has a generalized form 

^^^^^^ + dx = °' 

where k,{x) is a known function of the space coordinate, denoting effective "capacity" (e.g. 
porosity of the medium) the traditional Godunov value of net function q at the moment 
to + At 

qi = <li - ^ [Fi - Fi+i) , 
where Fj is the flux of q from the ith cell to the {i-l)th cell should be replaced by 

(it = qi T— [Fi - Fi+i) , 

where Axj is the size of the ith cell, while Ki is the cell-averaged value of k(x). 



2.5 Simulation Results 

Evolution of an accreting flow has been simulated according to the scheme described above 
for a number of sets of global parametres. 

Strong shocks have been found to evolve in the column on timescales of about 10~^ s. 
The shocks oscillate around their equilibrium positions with periods of about 10~^ s, the 
oscillations dumping time being about 10~^ s. A typical evolution of a shock and of a 
flow velocity profile is presented in Fig|21 

Modeled flow profiles in Figs. 01 and El demonstrate stable and strong shocks that 
decelerate and heat the accreting flow. At such shocks the ions are heated much more 
than the electrons as they contain most of the kinetic energy of the flow. However, as the 
hot ions move towards the surface in the downstream zone, they pass a substantial part 
of their energy to the electrons. In their turn, the electrons emit the energy as cyclotron 
and Bremsstrahlung photons and pass it to non-resonant photons in Compton collisions. 

In most of the modeled cases the compression ratio at the shocks slightly exceeds 4 
(the maximal value for non-relativistic single- fluid shocks) due to weak relativistic changes 
in the adiabatic index of the ions heated upto a few tens MeV. 

An important property of the model is transformation of a substantial part of the ram 
energy of the flow into cyclotron photons of the optically thick line. Fig|71demonstrates the 
dependence of the part on magnetic field strength. The pressure of the trapped cyclotron 
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photons noticeably affects the stopping of the accreting flow. The cyclotron cooling of 
the ions is significant at magnetic fields exceeding 5- 10^^ G; in this case, the particular 
accretion regime significantly depends upon the detailed structure of magnetic fields in 
the star's atmosphere about 10^ cm from the surface. At such heights the field may be 
significantly non-dipolar due to the presence of local high-multipole components. The 
non-uniform structure of the field makes the column transparent for optical and X-ray 
cyclotron photons emitted by ions and electrons. Accretion regimes for such an optically 
thin column where the cyclotron emission of the ions freely comes out are demonstrated 
in Fig|3 An account for cyclotron cooling of the ions at magnetic fields exceeding 5- 10^^ 
G leads to a significantly lower position of the shock. In this case the emission spectrum 
will most probably contain a prominent proton cyclotron line in the optical/UV range. 
Due to strong deceleration and efficient cooling of the flow in an optically thin part of the 
column only about half of the flow energy reaches the surface of the star (see Fig|HI). 

If magnetic field structures are regular at the heights of about 10^ cm from the surface, 
the proton cyclotron line may become optically thick. In this case the accretion regime 
is analogous to that in Fig|2l as the high-frequency collisionless relaxation recovers the 
isotropy of the ion distribution much faster than the transverse temperature changes. 

It should be also noted that the qualitative difference in flow profiles at low and high 
values of magnetic field (Fig 121) is due to the fact, that at low fields the electron tem- 
perature in the downstream region is lower and the gradient of locked cyclotron photons 
energy density increases. A large gradient makes the flow to decelerate more efficiently. 

3 Destruction of CNO Nuclei 
in an Accretion Column 

An important characteristic of an accreting flow is the chemical composition of the matter 
that reaches the surface of a NS. The characteristic is particularly important for the theory 
of X-ray bursts, as accreted CNO nuclei may catalyze thermonuclear burning of hydrogen 
on the surface of a NS (see Lewin, van Paradijs, and Taam 1993, Strohmayer and Bildsten 
2003, and references therein). 

Most of the bursts (Type I bursts) have been observed in low-mass X-ray binary 
systems where the NSs are usually thought to possess magnetic fields not exceeding 10^ 
G, while the model presented here treats much stronger fields usually attributed to high- 
mass X-ray binary systems. Nevertheless, it is instructive to study the dependence of the 
survival fraction on magnetic field and mass accretion rate in the scenario of collisionless 
matter stopping at an accretion shock wave described above. 

The nuclei spallation process leads to fast destruction of tens MeV/nucl range nuclei 
and thus, to a significant suppression of gamma-ray line emission from accreting objects 
(see e.g. Aharonian and Sunyaev 1984 and Bildsten et al. 1992). Bildsten et al. (1992) 
considered the destruction scenario for the case when the accretion flow is decelerated by 
Coulomb collisions in a dense atmosphere of a NS and concluded, that almost all the CNO 
nuclei would be destroyed before they reach the surface. The authors noted, however, that 
this conclusion may be not true in the case when the flow decelerates in a collisionless 
column above the atmosphere as it is in our model. 
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Indeed, in that case, the depth traversed by a nucleus as it decelerates down to energies 
of about 10 MeV/nucl can be much smaller than the depth for a pure Coulomb deceleration 
down to the same energies. Having modeled flow profiles in an accretion column, one is 
able to say how efficient the destruction of CNO nuclei in an accretion flow is and where 
they are destroyed. 

In order to test the destruction efficiency quantitatively the destruction probability 
for a carbon nucleus that accretes inside a modeled flow was calculated (the destruction 
crossections of nitrogen and oxygen are very similar to that of carbon, so they would be 
destroyed in the same way). 

The carbon nuclei destruction rate [with crossections given by Read and Viola (1984)] 
was numerically integrated within the flow to obtain the nuclei survival fractions as a 
function of the distance from the NS surface. 

The dependencies of survival fractions for a carbon nucleus on the distance from the 
NS surface are presented in FigE] for a set of magnetic fleld values. It follows that the 
survival fraction of C nuclei drops by an order of magnitude for magnetic flelds increasing 
from 5-10^^ G to 5-10^^ G. For magnetic flelds about 5-10^^ G a signiflcant fraction of 
the nuclei may reach the surface of the star. The survival fraction is not too sensitive to 
the accretion rate in the interval studied (from 10^^ g s~^ to 10^^ g s^^). The effect may 
prevent Type I bursts for NSs in HMXBs with magnetic flelds > 3 ■ 10^^ G, while the 
bursts are allowed for lower flelds, if the hydrodynamic accretion regime considered in the 
paper realizes. 

4 Summary 

A numerical model of a non-stationary one-dimensional accretion column over a polar 
cap of a magnetized neutron star has been constructed. The model describes a two-fluid 
flow of accreting plasma in a strong dipolar magnetic fleld. One of the main features of 
the model is an ability to treat consistently flow discontinuities (including shocks) with 
a modifled Godunov method that made it possible to study temporal evolution of shocks 
in the accreting flow. 

After several free-fall periods a quasi-stationary state of the column with a stable 
accretion shock is usually reached. At the shock the ion temperature jumps upto ~ 
10^^ K (this value weakly depends upon the magnetic fleld strength and accretion rate). 
Depending on the magnetic fleld strength the electron temperature reaches (3-5) x 10® K. 

A part of the kinetic energy of the flow is transformed into emission of a thick cyclotron 
line. The radiation pressure in the line signiflcantly affects the deceleration of the plasma 
flow. A substantial part of the kinetic energy of the flow is emitted into the optically 
thin part of the spectrum well before the flow reaches the bottom of the column. It is 
often assumed that the kinetic energy of the flow and the emission of the optically thin 
part of the column are transformed into blackbody radiation in the optically thick part 
of the star's atmosphere. As a signiflcant part of the flow energy is emitted far above 
the atmosphere the effective temperature of the polar cap T^ff is reduced and the usual 
relation T^jf oc M^/^ is not true anymore, because the part of flow energy emitted above 
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the surface is now a complicated function of M. 
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Appendix: Energy Exchange Rates 



The rates of energy exchange due to small-angle Coulomb collisions ifgi (and Hi^ 
—H^i) were adopted from the work of Langer and Rappoport (1982): 



-A, (17) 



where A is a Coulomb logarithm, while Tq = AirrlrrieC^ is a typical energy-loss 
timescale. 

2. To calculate the coohng rates Qcyc and CyCei the complete QED crossection of 
coUisional excitation of electron Landau levels in a strong magnetic field adopted 
from the work of Langer (1981) was integrated numerically. In order to ensure 
reasonable accuracy in a wide range of temperatures and magnetic fields, ten Landau 
levels were taken into account. 

3. To calculate CyCee the following approximation was adopted from the work of Langer 
and Rappoport (1982). 



Cycee= 2.04ron^S-'/y^x (18) 
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Tfl C ! / 
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4. To calculate the Bremsstrahlung cooling rate Bvee we numerically integrated the 
crossection of Haug (1975) and constructed the following approximation: 

Sree ~ 2.5410 • 10-^X^-^^^^^^e • 9{B, Te), (19) 

where g{B,Te) = (0.409 - 0.01935i2 - 0.002445^2) (A;BTe/10 keV)°-^^ - is a gaunt- 
factor adopted from the work of Langer and Rappoport (1982). 

5. To calculate the Bremsstrahlung coohng rate Bvei we numerically integrated the 

crossection of Berestetskii and Landau (1982) for high electron temperatures and 
adopted the approximation from the work of Langer and Rappoport for low electron 
temperatures. 

_/0.36aro(Te/T,^)"-Vn.Z2(?(5,Te), when Te < T^ 
"""^^^ ~ \ OMaro{TjT^y-'neniZ'g{B,Te), when T^ > T^ ^^"^ 

where T^=5-10^ K. 

6. Q relax denotes Coulomb relaxation of the accreting flow on the electrons of the dense 
and thin atmosphere of a NS. The simplest model of an isothermal atmosphere is 
used and Qreiax is defined as 

Qrelax ^ei T(itm)i (^-i) 

7i - 1 

where Vei is the frequency of Coulomb collisions. 

13 



7. Qcompt denotes cooling of the flow electrons in single Compton scatterings. If a 
scattering photon is not too hard ^1), the energy lost by an electron in a 

single scattering is AE = — ^ + ^hsl^_ As < >= SksT^, < El >= UkeT^, 

Qcompt = ne'n^ < arVreiAE > H{B,T^) = (22) 

T — T 

= Unen^aTcksT^kB— :^H{B,T^), 

where H{B,T^) = (1 + 0.0l65{hujB/kBT^Y-*^)/{l + 0.0825(/icub/A:bT^)2-48) is the 
gaunt-factor adopted from the paper of Arons, Klein, and Lea (1987), while and 
are local values of photon temperature and energy density 
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Figure 1: Accretion column geometry 
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Shock formation and evolution 
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Evolution of an accreting flow: a) shock height vs time 

b) evolution of velocity profile 
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Flow Profiles at time = 1.3-10-3 g 
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Figure 3: Flow profiles at different magnetic field values 
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Flow Profiles at time = 0.8-10"^ s 
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Figure 4: Flow profiles at different magnetic field values when the column is optically thin 
for ion cyclotron emission 
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Flow Profiles at time = 5.0- IQ-^ s 
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Figure 5: Flow profiles at different accretion rate values 
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Figure 6: Survival fraction of ^^C nuclei on their way downto the surface 
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Figure 7: Part of flow energy converted into cyclotron emission at different magnetic fleld 
values 
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Figure 8: Temporal evolution of the part of flow energy that reaches the surface 
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